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Abstract 

We study the asymptotic behavior of Multiple Meixner polynomials of first and second 
kind, respectively [7] . We use an algebraic function formulation for the solution of the equilib- 
rium problem with constrain to describe their zero distribution. Then analyzing the limiting 
behavior of the coefficients of the recurrence relations for Multiple Meixner polynomials we 
obtain the main term of their asymptotics. 
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1 Introduction 

We consider the asymptotics of polynomial sequences (P n ) defined by orthogonality relations 
with respect to a discrete measure [i (with finite moments) 

N 

H = ^2p(x)5 Xk , p(x k ) > 0, x k G R, iVeNU{+oo}, (1.1) 

fc=0 

which is a linear combination of Dirac measures at the points xo, . . . ,xjy. By N we denotes the 
set of all nonnegative integers. 

The asymptotic theory of orthogonal polynomials with respect to discrete measures is not so 
widely developed as for absolutely continues measures of orthogonality. It can be explained be- 
cause the distribution of zeros of the polynomials orthogonal with respect to a discrete measure 
has an extra constrain. Indeed (it follows from the interlacing property of orthogonal polynomi- 
als), between two neighboring mass points of the measure of orthogonality it can not be situated 
more than one zero of the orthogonal polynomials. This constrain plays an important role in the 
logarithmic potential description of the limiting measure of zeros distribution and of the main 
term of asymptotics. 
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Logarithmic potential 

V(z) = - j \n\z-t\ dv(t), 

for probability measure i / p„(t) (zero counting measure), equally distributed in the zeros of poly- 
nomial P n , is 

V VPn = --lnlPJ . 
n 

Thus, the potential of a limiting (when n — >oo) zero distribution measure A defines an exponent 
of the main term of the polynomial sequence (P n ) asymptotics. Such type of asymptotics 
are also called weak asymptotics. For the description of this limiting measure it is useful to 
formulate a minimization problem for the energy of the logarithmic potential. It is clear, that 
constrain of the zero counting measure (which we pointed out above) leads to a problem of 
energy minimization in the class of probabilistic measures constrained by the weak limit of the 
mass points counting measure. For discrete orthogonal polynomials this approach was suggested 
by Rakhmanov in [21] . This approach was later extended in [9] (for more details on the topic we 
refer to monographs |19j and [22] for extremal and equilibrium problem of logarithmic potential 
with application to the weak asymptotics of polynomials sequences as well as to monograph [H] 
for discrete orthogonality and constrained equilibrium). 

The content of this paper starts with the study of weak asymptotics for monic polynomial 
sequences defined by discrete measure (jl.ip . In particular, we focus our attention in a polynomial 
sequence orthogonal with respect to negative binomial distribution (Pascal distribution) on N, 
i.e 

Pk = -^V, (B > , < c < 1) on x k = k, k G N. (1.2) 
k\ 

This orthogonal polynomial sequence (denoted by M n (x;/3,c)) is called classical Meixner poly- 
nomials [15] . They satisfy the orthogonality conditions 

^M n (£;;/3, C )(-£0;^V = 0, j = 0, 1, . . . , n - 1, (1.3) 

fc=0 

where (—k)j = —k(l — k) ■ ■ ■ (j — k — 1), j G N, and (k)o = 1, is the Pochhammer symbol. 
However, our main goal is the asymptotic analysis of multiple Meixner polynomials, which con- 
stitute a generalization of the aforementioned Meixner polynomials. For these polynomials the 
orthogonality conditions are considered with respect to a collection of Meixner weights (jl.2p with 
different parameters (3 or c. The study of these polynomials (among others classical discrete 
multiple orthogonal polynomials) was initiated in [JJ. 

The structure of the paper is as follows. The next subsection 11.11 is devoted to the notion 
of multiple orthogonal polynomials (definition and some properties). In subsection 1 1 . 21 we give 
general information about recurrence relations for multiple orthogonal polynomials. The intro- 
duction finishes with two subsections devoted to some formal properties of the multiple Meixner 
orthogonal polynomials obtained in [7J. Recall that there are two kind of multiple Meixner 
orthogonal polynomials. The first kind corresponds (see subsection II. 3p to simultaneous orthog- 
onality conditions with respect to discrete weights formed by (jl.2p with various q, i = 1, . . . , r. 
The second kind polynomials (see subsection ll.4p appear when the collection of discrete weights 
is formed by (jl.2p with various i = 1, . . . ,r. The explicit expressions for the coefficients of 
the recurrence relations of the multiple Meixner orthogonal polynomials are the main outcome 
from these subsections 11.31 and 11.41 -used in the sequel. 

In section [2] we state the logarithmic potential equilibrium problems for the description 
of weak asymptotics of Meixner and multiple Meixner polynomials and the solutions of these 
problems by means of certain algebraic functions are found. The procedure for solving such type 
of equilibrium problems is inspired in [2] (see also earlier references therein). 
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Lastly, in section [3] we, starting from the coefficient of the recurrence relations, obtain the 
main term of asymptotics of multiple Meixner polynomials and then we check the connection 
of this term with the equilibrium problem from the previous section. Theorems 13.11 and 13.21 are 
the main results of the paper. 

Furthermore, the spectral curves for the multiple Meixner orthogonal polynomials (of the first 
and second kind) and their connection with recurrence relations and vector potential equilibrium 
with external field and constrain constitute the main outcome of the present paper. 

Finally, we should also mention paper [23] by Sorokin in which the weak asymptotics of 
multiple Meixner orthogonal polynomials of the second type were studied, but using another 
approach where the integral representation for the polynomials and saddle point method were 
applied. 



1.1 Multiple orthogonal polynomials (type II) 

Below we summarize some results needed for the sequel. 

Assume one has r positive measures /ui, . . . , fJ, r on R, where for each i = 1, 2, . . . , r 

supp(^j) = jx £ R Ve > 0, Hi (^(x — e, x + e)^ > j , 

denotes the support of Hi- By fij we denote the smallest interval that contains supp(//j). 

We introduce a multi-index n = (m, n2, . . . , n r ) £ N r and its length \n\ = n\ + n,2 + . . . + n r . 
Multiple orthogonal polynomials (of type II) can be defined as follows |19[ Chapter 4.3], [I] and 



Definition 1.1. A type II multiple orthogonal polynomial P? L , corresponding to the multi-index 
ft £ W , is a polynomial of degree < \n\ which satisfies the orthogonality conditions 

/ Pn{x)x k dfii(x) = 0, k = 0, 1, . . . ,n» - 1, i = l,2, ...,r. 

Here one has a linear system of \n\ homogeneous relations for the |n| + 1 unknown coefficients 
of Pf[. Notice that the solution which is unique up to a multiplicative factor and also that 
this polynomial solution has exactly degree \n\ (then the monic multiple orthogonal polynomial 
exists and will be unique). When this happens n is said to be a normal index for type II. 

Suppose that the above r positive measures on R are discrete, i.e. 

Hi = ^2pi,A i>k , Pi,k > 0, x^f. £ R, iVi£NU{+oo}, i = l,...,r, 
k=o 

with all the x% k = 0, . . . , JVj, different and this for each i. In this case we have that supp(/Uj) 
is the closure of {xj,fc}^ ano - that Oj is the smallest closed interval on R which contains 
{ x i,k}k=o~ ^ ne corresponding polynomials are then discrete multiple orthogonal polynomials 
[7]. The discrete measures in this paper will be supported on N or a subset, which is achieved 
by taking x^fe = k for i = 1, . . . ,r. The orthogonality conditions are then more conveniently 
expressed in terms of the polynomials (-x)j. 

Definition 1.2. A discrete multiple orthogonal polynomial of type II on the linear lattice, cor- 
responding to the multi-index n £ N r , is a polynomial P^ of degree < \n\ that satisfies the 
orthogonality conditions 

^2Pft(k)(-k)j p ijk = 0, j = 0, 1,... ,th - 1, i = l,2,...,r. 

k=0 



3 



In what follows we will focus on a systems of measures 

N 

Hi = S ^Pi,kO~x k , p iik >0,x k eR, TV <E NU {+00}, i = l,...,r, 

k=0 

where supp(/j,i) is the closure of {xk} k=0 and 0, = O, for each i = 1, ...,r, for which every 
multi-index is normal, namely an AT system (see Nikishin and Sorokin |19l p. 140] as well as 
Arvesii et al. [JJ p. 26]). An interesting property related to the location of the zeros of the 
discrete multiple orthogonal polynomials is given by the following theorem. 

Theorem 1.1. (see [7], p. 26]) Suppose we have an AT system of r positive discrete measures. 
Then every discrete multiple orthogonal polynomial Pf t of type II, corresponding to the multi- 
index ft with \ft\ < N + 1, has exactly \ft\ different zeros on ft. 

Observe that in an AT system every multi-index ft, with \ft\ < N + 1, is normal. We refer to 
p. 26] for the following two examples of Chebyshev systems: 

Example 1.1. [19, Example 3 on p. 138] The functions 

v(x)cf,v(x)xcf, . . . , v(x)x ni ~ l cf , 

v(x)cr,v(x)xc^, . . . , v(x)x nr_1 c^, 

with all the Ci > 0, i = 1, . . . , r, different and v a function which has no zeros on M + , form a 
Chebyshev system on R + for every ft = (rti, . . . , n r ) S N r . 

Example 1.2. The functions 

v(x)T(x + /3 1 ),v(x)xF(x + /3i), . . . , v{x)x Tll - 1 T(x + ft), 

v(x)T(x + j3 r ), v(x)xT(x + fa), v(x)x nr ~ l r(x + /3 r ), 

with fii > and fii — (3j ^ Z whenever i ^ j and v a function with no zeros on IR + , form a 
Chebyshev system on HL + for every ft = (m, . . . , n r ) E N r . // /3j — (3j ^ {0, 1, . . . , N — 1} then this 
still gives a Chebyshev system for every ft = {n\, . . . , n r ) for which rij < N + 1, i = 1, 2, . . . ,r. 

1.2 Recurrence relation 

Suppose that all the multi-indices are normal for the r measures p±, . . . , p r . Then there exists an 
interesting recurrence relation of order r+1 for the monic multiple orthogonal polynomials of type 
II with nearly diagonal multi-indices |27j . Here the nearly diagonal multi-index, corresponding 
to n, is given by 

s(n) = (k + 1, k + 1, . . . , k + 1, k, k, . . . , k), 

s v y s v y 

s times r—s times 

with n = kr + s, < s < r. If we write P n {x) = Pg^(x), then the following recurrence relation 
holds: 

r 

xP n (x) = P n+ i(x) + ^2 a n,jPn—j 

(x), (1.4) 

with initial conditions Pq = 1 and Pj = 0, j = — 1, —2, . . . , —r. 
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A more general form of this recurrence relation is given by 

r 

x Pft{ x ) = Pfl+ei ( x ) + bft t oPfi(x) + ^2 bfijPfi-tf. (x), 

3=1 

3-1 ^ 

where e*j is the ith standard unit vector in W and vj = &r—k- I n the case r = 2 the 

fc=o 

recurrence relation for the polynomials with nearly diagonal multi-indices (|1.4p gives the relations 
(n+1, n) — > (n, n) — > (n, n—1) —> (n— 1, n— 1) and (n+1, n+1) — > (n+1, n) — )• (n, n) — > (n, re— 1). 
The first relation follows from the general case (m + l,^) — > (711,712) — >■ (ni,n2 — 1) — > 
(n-i — 1, ri2 — 1) by setting ni = n2 = n. To obtain the second one we set n\ = n and n 2 = n + 1 
and interchange the measures [i\ and //2- 

1.3 Multiple Meixner polynomials (first kind) 

Let fit, . . . , fi r be defined as 

^i = J2^ LS ^ 0< Ci <l, i = l,...,r, (1.5) 

fc=0 

with all the parameters q different. The support of these measures is N and we have Ai = . . . = 

A r = R+. lfp$ N\{0}, we define the functions , i = 1, . . . , r, in C°° (K \ -/3 - 1, -/3 - 2, . . .}) 

as 

^ w = ( r rw l) r(xti) if x€ RX «-!,-».-»,... }u {-ft -/J -1.-/J-2,...}), 

[ if X G {-1,-2,-3,...}, 

(1.6) 

with simple poles at -/3 - 1, -f3 - 2, . . .. If j3 £ N \ {0} we define 

which are functions in C°°. These functions can be written as w^(k) = ((3)f,c^/k\, k £ N, 
i = 1, . . . , r. By Example 11.11 we know that the measures fj,i, . . . , fj, r form an AT system which 
guarantees that every multi-index n = (n\, ... ,n r ) is normal for these measures. 

The monic multiple Meixner polynomials of first kind, corresponding to the multi-index 
ft = (ni, . . . , n r ) and the parameters j3, c= (c\, . . . , c r ), is the unique monic polynomial M^' c of 
degree |n| which satisfies the orthogonality conditions 

Y t M§' S (k)(-k) j wf(k) = 0, j = 0,...,m-l, i = l,2,...,r. 

fc=0 

In [7] the following Rodrigues formula was found 

c k Y k fQ , IW(x + 1)^/1 \ T(p + \K\+x) 



In particular, for r = 2, some straightforward calculations yield 

ni / \ n% 



r(/3 + |n|)r(x + i)' 



( Cl ^l)-m( C2 -l)« 2 (13) j k \ {j - k )\ 



cTc2 2 (/3) ni+n2 ^ <^ (-n 1 ) k (-n 2 )j-k \ c 



j=0 k=0 
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where 

m=0 n=0 

is the first of Appell's hyper geometric functions of two variables [10]. From this explicit expres- 
sion can be found the coefficients of the recurrence relation. 

Theorem 1.2. (see [7J) Multiple Meixner polynomials (first kind) Mn[% 2 2 satisfy the following 
recurrence relations 

where coefficients are 

b ni ,n 2 = ni(2cti + 1) + n 2 {ai + a 2 + 1) + ai/3, 

c ni ,n 2 = (n\(a\ + ai) + n 2 {a\ + a 2 )) (ni + n 2 + /3 - 1), (1.7) 
d ni ,n 2 = (/3 + ni + n 2 - l)(/3 + ni + n 2 - 2)(ai + l)(ai - a 2 )aini. 
with ai = an j a2 = _£2_ . 

1.4 Multiple Meixner polynomials (second kind) 

Assume that the measures fj,±, . . . ,fj, 2 are given by 

+°° (o \ k 
fc=0 

with < c < 1 and all the /3j different. The support of these measures is N. In [7J was found 
that every multi-index n = (m, . . . ,n r ) is normal for the above system of measures fj,i, . . . , fj, r , 
whenever — /3j ^ Z for all i ^= j. 
Here 

-*(»)-{ ^(fl'V+i) R\«-i, -2, -s,...}u {-ft -,J-i, -/>-*...}>, 

[ if ace {-1,-2,-3,...}. 

(1.9) 

The monic multiple Meixner polynomials of second kind, corresponding to the multi-index n 
and the parameters (3 = (/?i, . . . ,/3 r ), > (/% — /3j ^ Z for all i ^ j) and < c < 1, is the 

unique monic polynomial M^' c of degree |n| that satisfies the orthogonality conditions 
Y J M§ c {k){-k) j w^{k) = Q, j = 0,...,m-l, i = l,2,...,r. 

k=0 

For this multiple orthogonal polynomial the Rodrigues formula is (see [7]) 



\n\ r 



c \ 


c — 




for 


r = 


c 




c — 




c 


Y 



r^+i) -pr/ r(A) r(/3 i + n i + x) - 

11 Vr(ft + x) r(A + n<) ; r(x + 1) ' 



/ \ ni+n 2 

M*$«°(x) = (^J (/3 2 )„ 2 (/3 1 )„ 1 F 1 1 : g 



ni+n 2 / : (-ni);(-n 2 ,/3i+ni); 

1:1;^ I c— 1 c— 1 



(ft) :-;(&); 



c ' c 



(r \ (r \ sr v (- ra i)fc(~ n 2)j-fc(/3i + (— ) 



= ft ^ k\{j-k)\{p 2 )^ k (fa)] 
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where 

p q k 

+00 n (aj)r+ s n (bj)r n («,•)« _ . 

3=1 3=1 3=1 a: ?/ s 

P ' S=0 Il(a J -)r+.n(^)r 11(77). 
3-1 3=1 3=1 

with a= (oi, . . . ,ap), 6 = (61,. . . c = (ci, . . . , c fc ), a = (ai, . . . ,aj), /? = (/5i , - - . ,/3 m ), 7 = 
(71, . . . , 7 n ) are the generalizations of the Kampe de Feriet's series [53] which are a generalization 
of the four Appell series in two variables. From this explicit expression can be found the 
coefficients of the recurrence relation. 

Theorem 1.3. (see [7]) Multiple Meixner polynomials (second kind) M^f^ satisfy the fol- 
lowing recurrence relations 

<'frW = + b nun2 M^f c (x) + c^M^L^) + dn^M^f^ix), 

where coefficients are 

b ni ,n 2 = m(2a + 1) +n 2 (a + 1) + afii, 

c ni ,n 2 = a(a + l)^nin 2 + m(ni + - 1) + n 2 (n 2 + /92 - 1)J, (1-10) 

dni,n 2 = a 2 (a + l)ni(ni + /3i - l)(n! + /3i - /3 2 ), 

2 Equilibrium relations and algebraic functions 

Inspired in |21] and [9] we introduce here the main notions of the potential theory approach to 
the weak asymptotics (starting from the discrete measure of orthogonality) of the polynomials 
P n (x), orthogonal with respect to a discrete orthogonality measure (jl.ip . Usually one starts 
with scaling of the problem, i.e. 

x — > t : x = x(t, n), 

(for example, x = tn) such that the zeros of the polynomial sequence P(t) := P n (x(t,n)), 
remain in a compact set t S K <s C while n — > 00. This scaling produces a dependence of the 
orthogonality weight on the parameter n (varying weight): 

Pn {t) = p{x(t, n)) =: exp{-nV{t) + 0(1)}. 

Accordingly, a dependence on n also occurs for its discrete support 

{tk,n} ■ x k = Xk(tk, n , n), keN, 

for which one can define a counting measure for mass points of the support 

T n (t) := - 5(t-t k n ), 

k 

and the weak limit of this measure 

r(t) = lim r n (i) . 

n 



( a:b;c; \ 

if;!;!: , ■<■■!> 
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Under some technical conditions (for details see [21] and [9]) the limiting measure 



up n (x(t,n)) ->• A(t), 
minimizes the following energy functional (with external field V)) 

m--= J (v»(t)+±v(t))d»(t), (2.i) 

among measures from the class 

H : \n\ = 1, n < t. (2.2) 

The condition of boundedness of the measure /x by measure r reflects an interesting peculiarity 
of the discrete orthogonality -due to the interlacing of masses of the discrete measure and zeros 
of the polynomials. 

The potential of the extremal measure A satisfies the equilibrium relations with the presence 
of the external field 

k, t G suppA, 

= K, t€ supp(r - A) n supp(A) := S, (2.3) 
^k, t G supp(r — A). 

Set £, where the combination of the potential and the external field, is equal to a constant is 
called equilibrium zone. A set where the equilibrium measure A touches its constrain, i.e. measure 
r, is called saturation zone. As it follows from equilibrium relations (|2.3p in the saturation zone 
the combination of the potential and the external field is strictly less than the equilibrium 
constant k. 

To describe the weak asymptotics of multiple orthogonal polynomials an important notion 
of vector potential equilibrium was introduced by Gonchar and Rakhmanov in [11] and [12] . 
Let A := {Ai, . . . , A r } be a collection of compact sets in C and let 

D = (<!,.,)'', i- 

be a real symmetric nonsingular positive definite matrix. An additional condition on D to be 
compatible with (Ai, . . . , A r ) is that dij > whenever Aj n Aj ^ 0. 
For a vector of measures 

M = (jJ,i,...,fJr), supply C Aj, j = l,...,r, 

the energy functional is defined as 

r r 
i=l j=l 

where I(fii,fij) is the mutual energy of two scalar measures 

I(fj,i,Hj) = [ [ In- — l -—dfi i (x)dii j {i). 
J A, Ja 3 \x ~ t\ 

The extremal vector measure A, minimizing the energy functional (|2.4p among all fl (with fixed 
masses of the coordinates) possesses the equilibrium properties for j = 1, . . . , r 

uf( x ) = ± di jV^ x )h^ ^supp^AJ, (2 . 5) 
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Here the vector U x = (J7f\ . . . , Uf) is called vector potential of the vector valued measure X with 
respect to the interaction matrix D. 

External fields V = (V±, . . . , V r ) such that Vj(t) is a function on t G Aj, j = 1, . . . , r, and 
constraining for fl measures as follows 

t = (ri,...,r r ) : fij ^ Tj , suppr,- C Ay, j = l,...,r, 

can be incorporated in the vector equilibrium problem in order to develop an approach for 
multiple orthogonality with respect to varying weights or (and) with respect to discrete weights 
(in a similar fashion that [2TJ and [9]). Thus, the energy functional (|2.1|) of vector measure with 
external fields becomes (see (|2.4|) ) 



ear) := m + ViW^i® 



3=1 

and extremal vector measure A, minimizing this energy functional possesses the equilibrium 
properties (see (|2.5p ) for j = 1, . . . , r 

x G supp Aj , 

x G supp(rj — Aj) n supp Aj = Ey, (2-6) 




x G supp(rj — Aj) 



A saddle point extremal problem for functional (|2.6|) is considered in some applications 
(when positions of the components of A are not known). In such situation one performs the 
minimization of the energy functional by measure fl and maximize it by vector A. As a result 
the support of the equilibrium measure possesses the following symmetry relations 

= °_El±m, „ nSjlJ=1 ,..., r , ,, 7) 

on + on_ 

where d/dn± denotes the normal derivatives on the respective contours. This property is called 
S-symmetry after Stahl, who introduced it in [25] and [26] . 

To conclude this short survey on the logarithmic potential approach to the weak asymptotics 
we recall one general class of multiple orthogonal polynomials introduced by Nikishin in [17] , 
Such systems are defined as follows. Let o~\, o~2 be two finite Borel measures with constant sign, 
whose supports suppoi, suppc^ are contained in non intersecting intervals Ai, A2, respectively, 
of the real line R. Set 

d(a 1) a 2 )(x) = j da2 ^ dai{x). 
J x-t 

This expression defines a new measure whose support coincides with that of o~\. In general, let 
a = (<7i, . . . , o>) be a system of finite Borel measures on the real line with constant sign and 
compact support. Let A& = [afe,6fe] denote the smallest interval which contains the support of 
<7fc. Assume that A^ n Afc +1 = 0, k = 1, . . . , m — 1. We say that S = (si, . . . , s r ) = Af(a) is the 
Nikishin system generated by a, when 

si = t7i, s 2 = (<ti,<J2), • • • , s r = (cri, (<7 2 , • • • ,o r )). (2.8) 

In what follows (without loss of generality) we assume r = 2. So if we introduce the following 
notation for the Markov function (Cauchy transform) 



s(z) :-- 



ds(x) 
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then Nikishin system of two measures is: s% = a*i, S2 = cr$ o\. These measures have the same 
support. The Nikishin system is a useful model for multiple orthogonal polynomials with respect 
to a collection of measures (1*1,1*2) with the same support (A := supp/^i = supp/^)- For such 
type of systems an important role plays the analytic continuation of the function 

~r- ^2+ - 1*2- a /r, n \ 

■r := :=r- , on A . (2.9) 

Mi+ - Mi- 

For the pair (f*i,[*2) forming a Nikishin system (j2.8|) . (f*i,f*2) := (si> s 2) we have 

J" = 02 € F(C\A 2 ). 

The limiting zero counting measure for multiple orthogonal polynomial P^, n := (n,n) with 
respect to a Nikishin system of measures (s\, S2) 

is described by equilibrium problem (|2.5|) with the matrix of interaction 

B := ( _\ ~\ ) , (2.10) 

and A = Ai, where |Ai| =2, suppAi € Ai and | A2 1 = 1, suppA2 G A2. Here the measure A2 is 
the limiting zero counting measure 



->■ A 2 , j = 1,2, 

for the functions of the second kind 

iff) := [ P ^ ds ^ X \ j = 1,2. 
J z - x 

Since the original contribution [18] by Nikishin there have been published numerous papers 
devoted to the asymptotics of multiple orthogonal polynomial with respect to a Nikishin system 
of measures (for most general results and latest surveys see [1] and [3]). 

In this section, starting from the fact that components of the vector of discrete Meixner 
weights have the same support we assume that these measures inherit properties of the Nikishin 
system of continuous measures with the same support (which is not developed for the discrete 
measures yet). This assumption allows us heuristically to state an ad hoc vector equilibrium 
problem with Nikishin matrix of interaction, with external field (because of unbounded set of 
orthogonality) and constrain (due to discrete weight). Accordingly, at this stage we look for the 
solutions of this equilibrium problem in a form of spectral curve (among the simplest candidates 
satisfying the necessary conditions following from the equilibrium problem). Thus, as a result 
of the first stage of our analysis we have heuristically got equations of the algebraic curves 
(|2.26p and (|2.39p and have related their logarithms (rigorously) with certain vector equilibrium 
problems (see Propositions 12.11 and 12. 2p . We highlight that at this point spectral curves and 
equilibrium problems are not formally related to the multiple Meixner orthogonal polynomials. 

In section [3] we derive the same spectral curves from the recurrence coefficients of the mul- 
tiple Meixner orthogonal polynomials. Consequently, it now rigorously relates Nikishin type 
equilibrium with constrains to Meixner discrete orthogonality. 
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2.1 Potential problem and its solution for classical Meixner polynomials 

It this subsection, aimed to illustrate our approach in connection with the above discussed 
techniques, we take (as example) the classical Meixner polynomials (| 1 .3[) . First, we state the 
equilibrium problem for the limiting zero counting measure. Then, we show our approach, which 
is based on an algebraic function formulation, by solving this equilibrium problem in a quite 
direct way. 

The scaling of variable x in (|1.2p 

x =: nt, 

stabilizes the zeros of the polynomials (i.e. when n — too, the zeros stay in a compact) 

P n (t) := M n {x), 

which one considers now as a polynomial of variable t. Indeed, it is orthogonal to {t k }^Zo with 
respect to the varying weight (depending on n) 

PnW- (tn) , c , 

concentrated in the points {^}fc 6 N. Taking the limit 

lim \pn(t)\ 1/n = c* = exp { -t In - 1 , 

ra->oo I c J 

we see that the external field for extremal problem (|2.1|) - (|2.2p and equilibrium problem (|2.3|) 
has the follownig form 

V(t) = In - Ret . (2.11) 



We also have 



-^5[t--\-^l(t)=-.T(t), t G M. + , (2.12) 

71 k=0 ^ 



where l(t) is the Lebesgue measure with unit density (dt) on the semiaxis R+. 
Thus, the limiting measure X(t) of zero distributions of polynomials P n (t) 



satisfy the extremal problem (|2.ip - (|2.2p as well as the conditions of equilibrium (|2.3p with 
external field (|2.1ip and constrain t|2. 12 j> . 

Let us show our approach for finding the explicit representation of the measure X(t). We 
look for the Cauchy transform of measure A (or Markov function) 

h(z) := \{z) = J , z e C\suppA , (2.13) 

as function h(z), which be analytic on some two sheeted Riemann surface 1Z := (C^°\ C^ 1 ^) . The 
possibility of analytic continuation of the function h(z) to another sheet of the Riemann surface 
through the compact S can be seen from the second relation in (|2.3p after it "complexification" 
and differentiation. Thus, the function h := (ho,hi) on 1Z has some natural condition for its 
values at oo. For the main branch ho, because of (|2.13p (and normalization of measure |A| = 1) 

h(z) :=" + ••• , ^^oo (0) , (2.14) 
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and for its continuation hi (taking into account derivative of the "complexification" of external 
field (gZHD ) 

1 1 



h(z) 



In 



+ 



■ oo 



(i) 



(2.15) 



The density of the measure A'(i) is proportional to the jump of imaginary part of the function 
h(z) along R+. Moreover, it is clear from electrostatic consideration of the problem (|2.ip - (j2.2j) 
that at the neighborhood of the left end point of suppA, i.e. to the right from the origin, this 
jump has to be a constant -here the measure A reaches its constrain by the Lebesgue measure 
l(t). Such behavior can be modelled by representing the function h{z) as logarithm of some other 
function, say E(z), which on this part of the support of the measure A (i.e. in the saturation 
zone) takes negative real values. Hence, we are looking for h{z) in the form 



h(z) = \nE(z), 



(2.16) 



where E(z) is a rational function on 1Z with one zero and one pole (at the points on the different 
sheets of Riemann surface 1Z with the same projection point 0) 



E(z) 



0, 

oo 



o(°), 
ow. 



(2.17) 



We have chosen the position of these points in order to provide holomorphicity of branches h 
-equation (|2. 16[) in C\M+. The conditions (|2.14p and ()2. 15[) give the following behavior for E(z) 
in the neighborhood of infinity on both sheets of 1Z 



E(z) 



1 + 
1 

c 



1 



1 



•oo(°), 
•oo«. 



(2.18) 



Observe that it is enough to have conditions (|2.17p and (|2.18p for determining E(z) and equation 
of its Riemann surface. Indeed, the Vieta relations give 



E 2 + 



1 ( 1 

z\c 



1 



1 



+ 1 



E + 



1 



0. 



From here it follows that Riemann surface of the function E(z) has branch points 



1- 1 
ei = — — 7= , e 2 = — . 
1 + V c ei 

Moreover, both branches of the function E(z) are negative on [0, ei] and positive on [e2,oo]. 
Therefore, the jump of the imaginary part of the function h{z) -see equation (|2.16p ~ is a constant 
on [0, ei] and equal zero on [e2,oo]. Thus, suppA = [0,e2], the interval [0, ei] is saturation zone, 
the interval [ei, e?\ = £ corresponds to the equilibrium zone, and [e2, oo] is a zone free from zeros 
of the scaled classical Meixner polynomials. Measure A obtained by this procedure satisfies the 
equilibrium condition (|2.3p . Furthermore, after integration one can obtain an explicit value of 
the equilibrium constant in (|2.3p 

, , l ~ c 
k=1 + In . 
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2.2 Potential problem and its solution for multiple Meixner polynomials 
(first kind) 

Assume r = 2, and n = (n,n). We rewrite weights (|1.6p as follows 

<(x) = ^±^e-^, , = 1,2. 
J x\ 

Here the supports of both discrete orthogonality measures coincide. Therefore, we pose an 
equilibrium problem for the weak asymptotics in a similar way to the Nikishin case. Indeed, we 
consider the equilibrium problem associated with the weights w^(x), and u{x)w^(x), where 



Br \ -, n till 

W%[X) _ xln^ _ 1 f e 



£2 

u(x) = "" lx "' J = e^" 1 ^" = — — / dt. 

v4{x) 2™ J r t-x 

The integral representation in the right hand side is due to the Cauchy integral formula with 
a contour of integration T encircling R + , i.e. it comes from infinity through the upper half 
plane along the positive semiaxis turning counterclockwise around zero in the left half plane and 
returning to infinity along the positive semiaxis through the lower half plane. In addition, the 
contour T serves as a boundary of holomorphicity region of an analog of the function (|2,9p . 

Now we do a scaling of the problem in order to compactify zeros of multiple Meixner poly- 
nomials (first kind). Let 

M^(x) = C n M^(xn) = x 2n + -... 

Recall that the function of the second kind associated with the multiple Meixner polynomials 
of first kind is given by 

r /3;c 



J^(*) = /^%-, 3 = 1,2 



where \Xj is the discrete orthogonality measure (HSJ). Let us denote R { J\z) = (nz). General 
feature of the Nikishin systems is that the functions of the second kind have ro-zeros on a 
compact of the complex plane outside of the support of the orthogonality measures. For the 
classical Nikishin system (see (|2.8p ) they accumulate along interval A2. For the case under 
consideration we expect that the function R$ {z) has n-zeros accumulated along an analytic arc 
7 in the complex plane. We assume that contour T contains this analytic arc 7. 

For the formulation of the equilibrium problem, we follow the Nikishin system model, in 
which the zeros accumulate along an open contour 7 C C (for the functions of the second kind) 
and in M + (for the multiple orthogonal polynomials). Since the mass points x^ are uniformly 
distributed over a subset of E, then the zero counting measure t/]^P;Sf x \ ~* -V where |A| =2, 
and A < I, being dl(x) = dx the Lebesgue measure, and supp(A) C R+. On the other hand, 

is equal to 1 with supp(/z) C 7. 1 

|A| = 2, \<l, supp(A) C 

H = 1) supp(^) =: 

Now we write the equilibrium problem for these measures (A, fj,) := A as (|2.6|) . where: 



. A fj,, and its magnitude is equal to 1 with supp(/z) C 7. Accordingly, 



(2.19) 

\fj,\ = 1, supp(^) =: 7 C T. 



The matrix of interaction D in the definition (|2.5p of the vector potential has the 
Nikishin's form (pTTO]) . 

Scaling of the weights wf and u brings external fields V\ := (— lnci)Rex on M.+ and 
V2 ■= f-ln^J Rex on T, respectively. 
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Thus, we have the equilibrium relations: 

Wi := 2V X -V 11 - (In c x ) Re x 
W 2 := 2V - V x - (in 2a) Rex 



< Kx, on supp(A) C R+, 

> K\, on M + n supp(/ \ A). 
= k 2 , on supp(/i) =: 7 C T, 

> k 2 , on r. 



(2.20) 



We also have to add to this system the ^-symmetry property fj2.7[) for the second plate F of our 
condenser (R + ,r): 

dw 2 dw 2 

on 7. (2.21) 



Now we look for a solution of the equilibrium problem (|2.19p . (|2.20p . (|2.2ip . In accordance 
with our approach (by differentiating) we pass from (|2.20p to the following relations for the real 
part of the Cauchy transforms of the measures A and /i, 

Re |2A — /2 + lncij = 0, on supp(A) n supp(/ \ A) =: S, 



Rc 



1 2/1 — A + lnc2 — lnci| = 0, on supp(/^) := 7, 



respectively. From these expressions one rewrites 

A = fi — A — In ci , on £ , 
/I + In C2 = A — ju + In c\ , on 7 . 

Notice that if we define on the three sheeted Riemann surface 



(2.22) 



U := (^o,^i,^ 2 



n 

7l 2 



C\E, 

c\{sUt} 

C\7, 



a function 



A, in ^0, 

iJ(^) := < ^-A-lnci, in Tlx, 
k -/2-lnc 2 , in 7^-2, 

then due to (|2.22p and (|2.2ip the branches of this function have analytic continuation one to 
another through the cuts joining the sheets (T^-cb T^x-> ^2) °f the Riemann surface 71. Observe 
that from (|2.19p follows that H(z) behaves near infinity as 

-1 



H(z) 




as z 
as z 
as z 



00 in 7^0) 
00 in 72-1, 
00 in 1Z 2 . 



In the same way like we did it for the classical Meixner polynomials we model the constrain 
condition A < I and dl{x) = dx on R+, representing the function H(z) as logarithm of a 
meromorphic function <f>{z) on TZ, i.e. 



H(z) = \ix(j)(z) 



G M(K). 



(2.23) 



To have a jump for the function H along R + lifted to sheet TZq we force the function 4>(z) to 
take values <j> = 00, and = at the points from the two sheets (TZq,TZx) which have on C 
projection point z = 

, as z — > in TZq , 
00, as z — > in 7£i. 



(2.24) 



14 




Thus, H €M\K \ {R^ U }) • From tlie behavior of H at infinity we conclude that 

as z — > oo in T^o, 
[z) = { cj" 1 (1 - z" 1 ) + ••• , asz^oo in fci, (2.25) 

as z — > oo in 7^2 • 

If £ and 7 are connected then the genus of the Riemann surface 1Z is equal to zero (we assume 
that it is true). The simplest meromorphic function <j) (i.e. see (|2.24p -one pole and one zero) 
that maps 1Z to the Riemann sphere and therefore the inverse function z((f>) is just a rational 
function. From (|2,25|) we know the poles z as a function of (j> and moreover we know its residues 
in the poles. Taking into account also (|2.24|) one gets that the algebraic function <p(z) is defined 
by the equation 

: z = -r-i ~ -x— r " -x— r- (2 - 26) 

0—1 C\<p — 1 C20 — 1 

This algebraic function has four branch points. In order to construct (for the condenser with 
two plates (K+,r), T C C \ R+) a solution of the equilibrium problem (|2.20p - (|2.2ip we need to 
have the case when: 

• 4> has two complex conjugated branch points e, e, and two real e%, e<i (0 < e± < 62); 

• the analytic arc 7 joining points (e,e), chosen by means of (j2.21[) . does not cross R+; 

• values (j>(() G R_, when ( G [0,ei] (0) C Kq. 

The first item can be checked by analyzing the critical values c\ , C2 when the branch points 
of the function <j) coincide -the transition of two conjugate branch points to the case of two 
real branch points. For this purpose we take consequently two times the discriminant of the 
curve (|2.26p -first time in the variable <j> and second time in the variable z. In result we get a 
polynomial of two variables c\ and C2 with rational coefficients. Factorization of this polynomial 
gives a factor of high degree (maximal degree in each variable Cj is 6 and maximal degree of the 
monomials is 9). It defines an algebraic curve T(c\,C2) = 0, which has genus equal 1. Thus, 
this curve is rather complicated. However, taking into account the symmetry of this curve with 
respect to (01,02), and changing variables 

ci + c 2 =: s , 
ci c 2 =: p, 

yields the algebraic curve 

f(s,p) := T{c 1 {s,p),c 2 (s,p)) =0, (2.27) 

which has genus equal 0. Therefore, curve T(s,p) admits a rational uniformization. Using 
Grobner basis technique we arrive to the following simple parametrization of this curve 

(t 2 + 8t + 18)(t + 2) 2 



2t 

(t + 3) 3 (t + 2) 3 

An analysis of this curve shows that c±, C2 G R when t G [—00, — 2] n [0, 00]. On Figure [Q (A) for 
t G [—4.73, —3] we see branches of the critical curve (|2.27|) which belong to the interesting for 
us quadrate (ci,C2) G {(0,1) x (0,1)}. If parameters (ci,C2) are inside of the domain bounded 
by these branches and axes {c\ = 0, C2 G (0, 0.5)} and {c2 = 0, c\ G (0, 0.5)} then function cf> has 
two complex conjugated branch points e,e, and two real ei,e2 (0 < e\ < €2)- We denote this 
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domain QM . On the critical lines the conjugate branch points e, e coincide on the real axes and 
when (ci, c 2 ) G {(0, 1) x (0, 1)} \ QM =: A, the function d> has four real branch points. 

The second item, i.e. condition on the parameters (01,02) implying 7 n = can be 
clarified using an equivalent (constructive) definition of the S-curve 7 (see (|2.21|) ). It is known 
(see [2], [25]) that S-property condition (|2.2ip means that the analytic curve 7 belongs to the 
union of the critical (orthogonal) trajectories of the quadratic differentials 



Re / (Hj - H k ) dz = 



(H j -H k ) 2 (dzf < 0. 



(2.28) 



Using uniformization (|2.23p . (|2.26p we perfom integration by parts 

(0-1) 2 



H(z) dz = z In 1 



z hid) + 



(4> - a 1 )(4> - a 2 )' 



l/cj, j = l,2. 



Thus condition ()2.28|) is equivalent to 



exp <^ / (Hj - H k )dz 



1 



<?i)(4>k - 0-2) 



D 2 ( 



a 1 )(4> j - a 2 ) 



1. 



(2.29) 



Now we obtain the critical values of the parameters such that curve 7 is crossing the origin. 
Substituting in (pT29|) j = 1, k = 2, z = and taking from (pT!H|) . (pT26j) value of d>(0) we have 



(&(0)-<7l)(&(0)-<7 2 ) 



and 



(02(0) -l) 2 

Therefore, taking into account that cr±, a 2 G K we arrive to 



(0) 



2<7icr 2 



CT 2 



(J\ + £72 



(C7l - 1)(<7 2 " 1) 



(CTI - CT 2 ) 5 



1 

4' 



or 



(<Tl + £T 2 - 2) 



(*1 



0"2j 



±1. 



^1 



0"2 



A branch of the critical curve corresponding to sign (— ) in the right hand side lies out of the 
admissible range for c\,c 2 (or a\, 02) ■ Finally, considering in the right hand side sign (+) we get 



(<ri +a 2 -2) 2 = 2(at -a 2 ) 2 . 



(2.30) 
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Figure 2: The sets 7 = supp(/i) and [0, e 2 ] = supp(A) as well as the branch points (circles) 
are depicted for the cases: (A) -on the left. (01,02) = (0.5,0.25) G M\ (B) -on the right. 
(ci,c 2 ) = (0.5,0.146445649) is critical. 



On Figure Q] (B) the branches of the critical curve (|2.30|) forms the boundary of the domain M. 
In this domain (01,02) G M the analytic arc 7 does not intersect (see Figure [2] (A)). On the 
boundary of this domain 7 is passing trough the origin (see Figure [2] (B)). 

The cut joining the logarithmic singularities of the function H starts at point on sheet TZq, 
goes along R + to branch point e\, passing to sheet 1Z\ and goes there back to point 0. Along 
this cut (for all (01,02) G {(0, 1) x (0, 1)}) values (j) are negative. Thus we have, indeed 

HeM(n \ {[0,ei](°)U[ ei ,0]«}) • 
Summarizing we have that for the parameters belonging to the domain 



( Cl ,c 2 ) G {(0,1) x (0,1)} : 



c 2 (l - yg) 
2c 2 - 1 - y/2 



< ci < 



c 2 (l + yg) 
2c 2 - 1 + V2 



(2.31) 



the density of measure A is determined by the limiting values at the points [0, e 2 ] of the imaginary 
part of In 4>(z)/ 'it or equivalently, 



dMx) .. 1 T rr , . . 
— - — = lim — 1m ti [x + ly) 

ax 7T 



x G [0, e 2 ] in TZq 



(2.32) 



The interval [0, e 2 ] is the support of measure A. Since the function (j) is negative on [0, ei] in 
TZq we have that the density of the measure A is equal to unity there. Therefore, the interval 
[0, ei] corresponds to the saturation zone and the interval S = [ei, e 2 ] is the equilibrium zone. 
The support of measure \x is the analytic arc 7. Thus, we have proved the following result. 

Proposition 2.1. The density of the absolute continues measure \, solving equilibrium problem 
()2.19p . ()2.20p . (|2.21|) for (ci,c 2 ) G N, is equal to (|2.32[) . where the function H is defined by 
(I2T23D and (I2T26I) . 

Remark 2.1. We show below in section \3.S\ that the logarithmic potential of the measure A 
defined in (|2.32p is equal to the exponent of the main term of the asymptotics of the multiple 
Meixner polynomials for all values of parameters (ci,c 2 ) G {(0,1) x (0,1)} (see Theorem \3.2\) . 
Correspondently (for all values of the parameters) the zero counting measure of the polynomials 
converges to A 



17 



and the interval between the origin and the smallest real branch point will be the saturation 
zone. However, for (01,02) 6 {(0,1) x (0,1)} \Af the measure A is not any more a solution 
of the equilibrium problem (|2. 19[) . (|2.20p . (|2.2ip with Nikishin's matrix of interaction. Without 
going into details, we just mention that measure A from (|2.32j) for (01,02) £ A can be presented 
as a solution of an equilibrium problem with Angeleso matrix of interaction (see }1 Jjj ). For 
{c\,C2) € QN\N it will be more sophisticated equilibrium problem for the condenser with three 
plates (see examples in J3]/, J^]/J. It is interesting to mention, that if we take multiple Laguerre 
polynomials considered in 115^ as the direct continuous analogue of the multiple Miexner 
polynomials we will see a complete correspondence of these three asymptotical regimes (M, A and 
QN\N) between the continuous and discrete case. 

2.3 Potential problem and its solution for multiple Meixner polynomials (sec- 
ond kind) 



Let r = 2 and < @2 < Pi < + 1 in the representation (|1.9p for weights of the discrete orthog- 
onality measures. Now writing the measures of orthogonality for multiple Meixner polynomials 
(second kind) in form of a Nikishin system, we consider the ratio of the weights 

w^(x) ^ u k 1 (fo -Jh + l) k 

This meromorphic function can be understood as the Cauchy transform of a discrete measure u 
with support belonging to R_. Thus orthogonality measures of multiple Meixner polynomials 
(second kind) form a Nikishin system generated by two discrete measures: fi\ (see (|1.8p ) with 
support in M + and u with support in R_. Since the mass points for both measures are uniformly 
distributed over a subset of M+ and R_, respectively, then the limiting zero counting measures 
in the equilibrium problem will be constrained by dl{x) = dx the Lebesgue measure. 
Doing a scaling 

M^ 2 ' ,c {x) : = M n^ r,c ( nx ) , R^ uP2 ' c) ' j (z) = R^ 2 ' c ^ j (nz) j = 1, 2 , 

we get the varying weight of orthogonality which defines an external field on M + . 

Thus, we expect that the zero counting measures associated with the scaled multiple orthog- 
onal polynomial and the functions of the second kind weakly converge 



M 2n lt n 



to measures (A,//) := A 



|A| = 2, \<l, supp(A)cM+, 

M = 1; A 4 ^ h SU PP(^) C R_. 



Now we write the equilibrium problem for these measures A = (A,/i) as (j2.6|) . where: 

• Matrix of interaction D in the definition ()2.5p of the vector potential Uj has the Nikishin's 
form (l2TTUj) . 

• Scaling of the weights u^ 1 brings external fields V\ := (— lnc) Rex on M.+. 
Thus, we have the equilibrium relations: 



H/i := 2k — T 7 ^ — in c Re x 



VF 2 := 27^ - P ; 



> «i, on n supp(/ \ A). 
< on supp(^) C E_, 

> ^2, on R_ n supp(/ \ /i). 



(2.34) 
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Now we look for a solution of the equilibrium problem (|2.33|) . (|2.34p . In accordance with 
our approach (by differentiating) we pass from (|2.34p to the following relations for the Cauchy 
transforms of the measures A and n, 



A = ju — A — In c, on supp(A) n supp(7 \ A) =: T,\, 

ju = A — ju, on supp(/u) Pi supp(Z \ fi) =: T,^. 

We point out that if we define on the three sheeted Riemann surface 



(2.35) 



K := (Ko,K u K 2 ) : 



n 
iz 2 



c\s A , 

c\{s A us^} , 



a function 



in 1Z , 
in TZi, 
in TZ 2 , 



H(z) := I fi- A - lnc, 

then due to ()2.35j) the branches of this function have analytic continuation one to another throw 
the cuts joining the sheets (JZo,1Zi,TZ 2 ) of the Riemann surface 1Z. Observe that from (|2,33p 
follows that H(z) behaves near infinity as 

r 2 



H(z) 



+ 



1 



In c h 

z 



as z — > 00 in TZq, 
as z — > 00 in 1Z\ and 1Z 2 ■ 



In the same way like we did it in the previous subsections we model the constrain condition 
(A < I on and [i < I on R_, d/(x) = (ix) representing the function -ff(z) as logarithm of a 
meromorphic function 0(z) on 1Z, i.e. 



= In <f>£M(n). 



(2.36) 



To have a jump for the function H along M + lifted to sheet TZq and along M_ lifted to sheet TZi 
we force the function (j)(z) to have pole (f> = 00 of order 2 at the point z = on the sheet 1Z±, 
and zeros <j) = of order 1 at the points z = on the two sheets T^o an d 7^-2 



0. 

00 s 
0, 



as z — > in 7^o, 
as z — > in 72.1, 
as z — > in 1Z 2 . 



(2.37) 



Hence, H £ M (k \ {r^ U IJ mL 2) |) • From the behavior of at infinity we conclude 

, as z — > 00 in TZq, 

+ • • • , as z — > 00 in 72. 1 and 7Z 2 . 

We write expansions of branches (j) up 



2 

1 + - + 

z 



1 



1 



(2.38) 



Now we obtain an equation for the algebraic function 
to 0(l/z 3 ) 





= 1 + 


2 

- + 

z 


A 
~- 




= (7 — 


a 

- + 

z 


B 
~ 




= a — 


a 

- + 

z 


D 

T 1 
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Since is a rational function on the Riemann surface 1Z and all its zeros and poles are shown 
in ()2.37p . we consider the Vieta relations for the product and for the pairwise products of these 
expansions and obtain the following linear system for unknown A and B + D: 

[B + D)+3a + aA =0 _ Aa - 1 _ a 2 + 2a 

[l + a){B + D)-Aa + a 2 + 2aA) =0 ^ ' <j — 1 + : ~ a - 1 ' 

Hence, we obtain 

03 _ (2* 2 -2s-l)* + (* + !)» , + g (x-2) g + 2 (3; + l) _ ^ = Q 

This algebraic function has four branch points, namely infinity point oo and three real points 
e_,ei,e2 (e_ < < ei < e2). The cut joining the logarithmic singularities of the function H 
starts at point on sheet 1Z$ goes along IR + to the branch point e±, passing to sheet IZi goes 
there back via point to the branch point e_ and then again passing to sheet 7^2 goes there 
back to point 0. Thus we have, indeed 

H e M (n \ { [0, ei ](°) (J [ei,e_]« [j [e_,0p>}) . 

Recall that the density of the measure A is determined by the limiting values at the points [0, €2] 
of the imaginary part of the hi cj)(z) /n or equivalently, 

d\(x) 1 

— — = lim - Im H(x + iy) , x£ [0, e 2 ] in ^ . (2.40) 

dx y->-0 7r 

Since function 4> is negative on [0, ei] in 7?.o we have that density of the measure A is equal 
unity there. Therefore, [0, ei] corresponds to the saturation zone and E = [ei, is the equi- 
librium zone. Thus, we have proved the following result. 

Proposition 2.2. The density of the absolute continues measure X, solving equilibrium problem 
(|2T55]> . ([OlD is equal to ([2T3D]) . w/iere i/ie function H is defined by (gJSD and Q739]) . 

3 nth-root asymptotic from recurrence relations and equilib- 
rium problem 

In this section, starting from the coefficient of the recurrence relations, we obtain the main term 
of asymptotics of multiple Meixner polynomials and then we check the connection of this term 
with the equilibrium problem studied in the previous section. The main steps of our approach 
here are the following. 

1) We begin by forming a transition matrix for our recursion. We consider vectors 
and V 3 with coordinates 

M ni , n2 (x) 



(m + n 2 )! ' 

taken for the following values of the multi- index (711,712): 

[(n - 1, n - 1); (n, n - 1); (n, n)] ->■ V 1 , 

[(n,n-l); (n,n); (n + 1, n)} ->• K 2 , 
[(n,n); (n+l,n); (n + 1, n + 1)] -)■ K 3 . 



(3.1) 
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We define transition matrices An \ j = 1, 2, by means of relations: 

To obtain the coefficients of these transition matrices we need two types of recurrence relations 
for M nin2 , which connects the following indices 



/) (n,n) and (n + l,n), (n, n), (n, n — 1), (n — l,n — 1), 

JJ) (n + l,n) and (n + 1, n + 1), (n + l,n), (ra, n), (n, n— 1). 



(3.2) 



The first type of the recursions we already have. They follow from (jl.7p and (jl.lOp if we set 
(m, 712) = (n, n). To get the recursions of the second type we substitute (m, n.2) = (n, re + 1) in 
(|1.7p and (jl.lOp and then swap ai with a 2 and /?i with /?2, respectively. 

2) Analysis of recursions for polynomials (|3.ip with indices (|3.2p (these recursions follows 
from (|1.7p and (|1.1U|) ) shows that if we change the variable there as 

x 

t := - , 

n 

then we get the recurrence relations with coefficients having limits when re — > 00. It gives us a 
regime 



n 



00, 



(3.3) 



t g K m C, 



for investigation of asymptotics starting from the recurrence relations. Thus we form a limiting 
transition matrix 

A(t) = lim A® AW t G K m C, (3.4) 

n— >oo 

and find eigenvalues of this matrix. We fix the following order for these eigenvalues 

\Li{t)\ > \L 2 {t)\ > |L 3 (t)|, tGC. 

3) Then, nth. root asymptotics for polynomials (|3.1|) can be obtained by means of the following 
lemma. 

Lemma 3.1. There holds the following asymptotic formula (when n — > 00) 



I In 

n 



M n , n (x) 



(2n)! 



- / In 

n 



n 



dh + o(l) 



(3.5) 



uniformly for — E K <£ CI, where Q is a domain containing 00 -point and bounded by a curve 
T:= {t:L!(i) = L2(i)}. 

The proof of this Lemma follows from Poincare's theorem on ratio asymptotics of solutions 
of the recurrence relations (see [20]) and from the theorem of Kuijlaars and Van Assche on the 
nth root asymptotics of the solutions of the recurrence relations with varying coefficients (see 
[T3] as well as [S] and [B]). 

We note that function L(t) is an algebraic function since it is a root of the polynomials in L 
with rational coefficients in t (which is the characteristic polynomial of the matrix A(t)). This 
observation allows us to use uniformization of the algebraic curve L(t) 



L = F 1 (s) 

x 

n 



where s is the parameter of uniformization, 



F2(s) 
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to evaluate explicitly the integral in the right hand side of (|3.5p . 

4) Finally we check that the obtained explicit expression for the right hand side of (|3.5p is the 
logarithmic potential of measure A, which is a solution of the corresponding vector equilibrium 
problem from the previous section 











f In 


Mi) 


n J 







dn = T X (t) 



For this purpose we differentiate 



d_ 

dt 



In L 



n 



dh 



-■ H(t) , 



and we verify that 



tf(i) = ln0 o (*) 



d\(0 



where <j>{t) is the algebraic function from the previous section, which gives an explicit solution 
of the equilibrium problem. 



3.1 Asymptotics for multiple Meixner polynomials of the second kind 

We start with the polynomials of the second kind since the calculations for obtaining the asymp- 
totics look simplest. Indeed, the following result is valid. 

Theorem 3.1. The main term of asymptotics for multiple Meixner orthogonal polynomials of 
the second kind has the form 



n 



{2n)\ 



n 



- +o(l), 
xJ 



and the convergence is uniformly for — E K <e C, where at the right hand side stands the 
logarithmic potential of the equilibrium measure X for the equilibrium problem f|2.33[> . ([2.341) . 
Furthermore, the Cauchy transform of measure X has the form 

X(t) = H(t) = In Mt) , teC\[0,e 2 ], 

where 4>o(t) = 1 + j, t — > oo is branch (|2.38[) of the algebraic function (p, defined by equation 
(|2.39|) . and e 2 is the maximal positive branch points of the function <f>. 

Proof. Following the step 1) of our approach (see description at the beginning of this section) 
we compute the transition matrices An and An \ and we consider them in the limit regime 



(|3.3p . Doing it we note that the obtained limits coincide. This circumstance allows us to use 
the matrix 

A^(t) = lim A W , t e K m C , 

n— >oo 

instead of the transition matrix A(t) from (j3.4p . We just have to multiply by 2 the right hand 
side of formula (13.50. Hence, we have 



/ u 



A^\t) 



1 



■fa(o + l) 

1 



|a 2 (a+l) \ 
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The characteristic polynomial P(L,t) of the matrix A^(t) has the form: 

P = L 3 - - L 2 t + - L 2 a + -a 2 L + -aL + - a 3 + - a 2 . 
2 2 4 4 8 8 

The equation P(L, t) = defines an algebraic function L(t) of genus zero, therefore this function 
allows a rational uniformization. Moreover, taking into account the linear appearance of t in 
P(L, t) we immediately obtain this uniformization if we take L as the uniformization parameter: 

L = L, 

P(L, t) + \L 2 t _ 1 (2L + b + b 2 ) (2L-b + b 2 ) (2L - 1 + b 2 ) 



f 



±L 2 4 L 2 



where a is changed (a := ft 2 — 1). From here, since t = we have 

= 4L2x 

n (2L + 6 + 6 2 ) (2L - 6 + ft 2 ) (2L - 1 + b 2 ) ' 1 j 

Now we can integrate by parts the right hand side of (|3.5p . i.e. 

2 ^ ln(L) dn^) = 2n ln(L) + x ln(2L + b + b 2 ) + x ln(2L - 6 + b 2 ) - 2x ln(2L - 1 + b 2 ) . 

Then, we substitute x(n) from (|3.6p into the obtained expression and multiplying it by 1/n, as 
a result we obtain the right hand side of (j3.5[) 

2 
n 



ln(L)dn = 21n(L)- 



4 L 2 V V (2L-1 + 6 2 ) 2 

It remains to check, that the real part of F{L) is indeed the logarithmic potential of the equi- 
librium measure A of the problem (|2.33p . (|2.34p . It will be true if the derivative of F(L) with 
respect to t coincides with the function H(t) = lncj)(t), where (j>(t) is the algebraic function 
(|2.39p . To differentiate F we use the uniformization variable L. Thus, 

dt 1 { " fo(L) \ (2L- 1 + 62)2 y 
Finally, substituting 

, _ fl _ (2£ + fr + fr 2 )(2L-fr + fr 2 ) 
0:_e ~ (2L-1 + 6 2 ) 2 5 

_ 1 (2L + 6 + 6 2 ) (2L - 6 + 6 2 ) (2L - 1 + b 2 ) 2 

z : - * - 4 L 2 ; 

6 := = , a = c _1 , 

we get an equation for the function $>(z), i.e. 

z 2 (f) 3 - (j) 2 (2az 2 + z 2 - 2za + 2z - a + l) + zcr(zcr + 2z - 2a + 2)4> - a 2 z 2 = , 
which coincide with (|2,39p . Theorem is completely proved. □ 



23 



3.2 Asymptotics for multiple Meixner polynomials of the first kind 

Finally, we consider the Meixner multiple orthogonal polynomials of the first kind. 

Theorem 3.2. The main term of asymptotics for the multiple Meixner orthogonal polynomials 
of first kind has the form 



I In 

n 



Af£vf(x) 



(2n) 



n 



V" - +o(l), 



being the convergence uniformly for ^ G K C, where at the right hand side stands the loga- 
rithmic potential of a measure A defined by its Cauchy transform which has the form 

A(t) = H(t) = In , teC\[0,e 2 ], 



where 4>o(t) = 1 + 4, t — > oo is branch (|2.38p of the algebraic function <fi, defined by equation 
(|2.26p . and e 2 is the maximal positive branch points of the function <f>. 

For parameters (c l5 c 2 ) G A/" (see (|2.3ip ) the measure A is the solution of the equilibrium 
problem (12+91) . (f2T20l) . ([231]) . 

Proof. Firstly, we compute the transition matrices A^ and ^4^ 2 \ and we consider them in the 
limit regime (|3.3|) 



A&'>(t) = lim A$ , t 



n 



For the case of the Meixner multiple orthogonal polynomials of first kind these limiting matrices 
are different. We have 



/ it 



fai 



\a 2 



1 — \a\{a\ + 1) — 502(02 + 1) 



and 



A (2) 



/ ^ 



V 



1 





— 2°i 





1 



■|ai(oi + l)(ai - a 2 ) \ 


/ 



1 -\a\(a\ + 1) — |a 2 (a 2 + 1) 



a 2 (a 2 + l)(a 2 - ai) \ 












/ 



Then we compute the characteristic polynomial P(L, t) of the transition matrix -product A^A^ (t). 
Like in theorem 13.11 the equation P(L,t) = defines an algebraic function L(t) of genus zero. 
Rational uniformization of this algebraic curve has the following form 



L 



(s - a 2 (a 2 + l)(ai - a 2 )) (s + ai(ai + l)(ai - a 2 )) 
(ai(a 1 + l) + (a 2 (a 2 + l)) 2 : 



_ ((s + af (ai + 1) + a\(a 2 + 1)) (s + oi(ai + l) 2 + a 2 (a 2 + l) 2 ) (2s + (a 2 + ai + l)(ai - a 2 ) 2 )) 
((ai(oi + 1) + a 2 (a 2 + 1)) (s + a x (ai + l)(a x - a 2 )) (s - a 2 (a 2 + l)(ai - a 2 ))) 

with s as a parameter of the uniformization. From here, since n = tx, we integrate by parts the 
right hand side of (|3.5p 



J ln(L) (in = n ln(L) — In 



s + a\ + af + a\ + a| 
s + ai + 2a 2 + af + a 2 + 2a 2 + a\ 



x . 
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Hence, multiplying it by — we obtain the right hand side of (|3.5 
i j ln(L) dn = ln(L(s)) + In ( 



s + ai + 2aj + af + a 2 + 2a^ + a^ 

It remains to check, that the real part of F is indeed the logarithmic potential of the equilibrium 
measure A of the problem (|2,19p . (|2,20p . (|2.2ip . Thus, we differentiate F using the uniformization 
variable s. Accordingly, we have 

i7 := ln((A) := — F = = - In ± 5 ^ ^ =: m(r(s)). 

v; eft £t V s + «i + 2 Oi + af + a 2 + 2al + al) Ky " 

From here we get a rational uniformization of the curve 

4>(z) : 4> = r(s) , z = t(s), 

which bring us an equation for the function (f)(z): 

a2diz(f) 3 + (— 3a2a\z — za\ + a\ — zai + a,2)4> 2 + (—2 + z — 0,2 — d\ + 1za,2 + 2zai + "&a2a\z)§ 

-z{a 2 + l)(ai + 1) = 0, 

or equivalently 

2 ai + 1 Q2 + 1 

11 1 1 ' 



<p — 1 ai</> — 1 — a\ ai<\> — a 2 — 1 

which after substitution a\ = and a 2 = coincides with (12.261) . It remains to apply 

1 — c\ 1 — c 2 

Proposition 12.11 The theorem is proved. □ 

4 Conclusions and future directions 

In this paper we have investigated the nth root asymptotic behavior of the multiple Meixner 
polynomials of the first and second kind, respectively (see theorems I3.HI3T2]) , We have devel- 
oped an approach based on an algebraic function formulation in connection with some available 
techniques from logarithmic potential theory. Indeed, for each polynomial sequence we have 
posed the corresponding equilibrium problem for the limiting zero counting measure and then 
we solved it in a quite direct way. A detailed analysis of the asymptotic distribution of zeros 
of multiple Meixner polynomials was given, i.e. we looked at different regions in the complex 
plane for identifying the support of the extremal measure, the saturation zone, the equilibrium 
zone, and the free zone. 

Finally, other systems of discrete multiple orthogonal polynomials (for r = 2) might be 
asymptotically analyzed using the approach presented here (up to some adaptations) to pro- 
duce similar outcomes. 

Acknowledgments. The authors are grateful to their colleagues Van Assche, Sorokin and 
Tulyakov for useful discussions. 
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